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Abstract Most of our knowledge of the physical processes in distant plasmas is 
obtained through measurement of the radiation they produce. Here we provide an 
overview of the main collisional and radiative processes and examples of diagnostics 
relevant to the microphysical processes in the plasma. Many analyses assume a 
time-steady plasma with ion populations in equilibrium with the local temperature 
and Maxwellian distributions of particle velocities, but these assumptions are easily 
violated in many cases. We consider these departures from equilibrium and possible 
diagnostics in detail. 
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1 Introduction 

Radiation is often the dominant cooling mechanism for optically thin astrophys- 
ical plasmas, which means that it determines the energy budget. It also provides 
most of the diagnostics for plasma parameters such as density, temperature and 
composition. It is therefore necessary to understand the dominant collisional and 
radiative processes in the plasma in order to answer astrophysical questions about 
the heating or energy dissipation in the plasma. In most cases, the radiation arises 
from collisions between electrons and ions, but interactions of electrons with a 
magnetic field or radiation field can also be important. 

The subsections of the introduction briefly summarize the processes that dom- 
inate in most astrophysical settings, including the wavelength ranges where they 
are observed and their identifying signatures. In this section we emphasize ra- 
diative signatures relevant to microphysical plasma processes, such as differences 
between electron and ion temperatures, turbulence, and non-Maxwellian velocity 
distributions. 
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The most detailed diagnostics for the physical parameters of plasmas and the 
microphysical processes taking place are generally based upon atomic and molec- 
ular lines and continua. In Section [2] we turn to a discussion of radiative processes 
and the diagnostics that are available. In Section |3] we present the theory of line 
formation in the coronal approximation and describe the dominant collisional and 
radiative processes. Section |3] comprises a discussion of the factors that influ- 
ence the charge state, including the key ionization and recombination processes, 
the charge state in temperature equilibrium and the circumstances under which 
the charge state can become decoupled from the local temperature. The micro- 
physics that arise when the electron distributions exhibit strong departures from 
Maxwellian are introduced in SectonO we review the kinetic equations that de- 
scribe the evolution of the distribution function and the different formalisms that 
have been adopted for handling collisions. In addition we address the consequences 
for the heat flux in terms of saturation and de-localization, and for the excitation 
and ionization rate coefficients which affect the ionization state and, in turn, the 
radiative losses. The optically-thin radiative loss function itself is the subject of 
Section [6] together with its dependence on the ionization state and the electron 
distribution. In SectionQwe return to a detailed review of the observational signa- 
tures and diagnostics that provide evidence for the importance of non-equilibrium 
ionization and non-Maxwellian electron distributions in the solar atmosphere. Fi- 
nally, we present a summary of our review and look to the future in Section [S] 



1.1 Bremsstrahlung 

Bremsstrahlung is continuum radiation produced by an electron when it is accel- 
erated in the electric held of an ion. The spectral shape is Pv ~ e xp(—hv/kT). The 



mechanism is well understood from basic electromagnetic theory ( Rvbicki &: LightmanI 
Il979h . but re lativistic corrections are needed for very high temperatures and pho- 
ton energies ( Nozawa et al.|[2009l) . Bremsstrahlung dominates the X-ray continua 



of many astrophysical sources, though the continua due to radiative recombina- 
tion and 2-photon processes should not be ignored (R avmond fc S mith 1 97711 . and 
there could be a contribution from synchrotron emission in young SNRs. 

Bremsstrahlung emission in the X-rays generally arises from thermal plasmas, 
but b remsstrahlung is al so seen from beams of non-thermal electrons in solar 
flares ( Kontar et al.fioill) . Bremsstrahlung emission is also referred to as free- free 



emission, particularly when observed at longer wavelengths. For example, free-free 
emission is observed from planetary nebulae and H II regions in the radio, and it is 
especially valuable as a measure of the ionizing flux from the central star, because 
it is unaffected by reddening. 

The signature of bremsstrahlung emission is a smooth continuum with an 
exponential cutoff at ~ kT. For normal astrophysical abundances it will be 
accompanied (and energetically dominated) by spectral line emission unless the 
temperature is so high that the abundant elements are ionized to their bare nuclei. 
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1.2 Synchrotron and Cyclotron Emission 

The emission from electrons Ryrati ng in a magnetic field can be accurately pre- 
dicted from electromagnetic theory ( Rvbicki 8z Lightman|[l979l) . Relativistic elec- 



trons dominate the radio and X-ray synchrotron emission from SNRs, the Galactic 
Halo, AGN and jets from X-ray binaries. Non-relativistic cyclotron emission can 
be important in the solar corona and in accreting magnetic white dwarfs. 

Synchrotron emission dominates the radio emission of supernova remnants, and 
in the fast shocks in young SNRs it produces narrow filaments of X-ray emission. 
The sharpness of the X-ray filaments is used to derive li mits on the diffusion 
coefficient for energetic particles in the acceleration region ( Long et al.|[2003i) and 



show that the magnetic field is amplified well beyond the values expected for 
compression in the shock (jVink fc LamindliooS iBamba et al.ll2005h . 



Cyclotron and synchrotron emission are highly polarized, but turbulence ran- 
domizes the field directions and Faraday rotation can change the polarizat i on di- 
rection and depolarize the emission from an extended region. iBvkov et al.1 (|2009|) 
demonstrate how t urbulence will affect the X-ray polarization on small scales, and 
iDickel et all (|l99ll) have shown that the radio polarization indicates radial, rather 



than tangential magnetic fields near the edge of Tycho's SNR. Polarization maps 
in the radio provide a u nique method for observing the turbulent structure of the 



m tne raaio provide a u nique metnod tor observing 
galactic magnetic field (jHaverkorn &z Heesenllioli ). 



Synchrotron emission dominates the radio and X-ray spectra of pulsar wind 
nebulae (PWNe), jets from AGN and gamma-ray bursts. It is straightforward to 
determine the power law slope of the emitting electrons from the slope of the spec- 
trum. The ambiguity between magnetic field strength and the number of emitting 
electrons can sometimes be resolved based on spectral breaks due to optical depth 
or synchrotron cooling. 

The emission and absorption occur between q uantized Landau levels in the 
solar corona at radio wavelengths ( Dulk et al.lll97d) . in magnetic cataclysmic vari- 



ables in the optical, and in accreting neutron stars in the X-ray. The emission at 
harmonics of the cyclotron frequency can be used to determine the magnetic field 
strength. The lowest harmonics often are optically thick and the higher ones opti- 
cally thin. At the transi tion, the radiation can be strongly polarized. For example, 
iBrosius fc Whitd (|2006l) used radio measurements above the solar limb to obtain 



the magnetic field strength above a sunspot. 

The signatures of synchrotron emission from relativistic electrons are a power 
law spectrum and a substantial polarization fraction. Gyro emission from non- 
relativistic thermal electrons typically shows a spectral peak corresponding to a 
modest harmonic of the cyclotron frequency, with substantial polarization. 



1.3 Compton and Inverse Compton Heating and Cooling 

The interaction between a photon and an electron can transfer energy either way. 
As for bremsstrahlung and synchrotron emission, the physical process is well un- 
derstood (Rybicki & Lightman 1979). Hot plasma above an accretion disk will 
experience Compton heating by hard X-rays from the central source and Comp- 
ton cooling by softer photons from the disk. Energetic electrons can interact with 
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synchrotron photons produced by the same electron population (synchrotron self 
compton emission). 

In supernova remants, the energetic electrons can produce TeV gamma rays 
by inverse Compton interaction with the cosmic microwave background (CMB) 
or with locally enhanced IR or optical radiation. It is currently debated whether 
the gamma ray emission observed from several SNRs arises from inverse Compton 
emission by energetic electrons or from decay of pions produced by interactions 
between cosmic ray protons and dense ambient plasma. Consideration of the lower 
energy gamma rays observed by FERMI can help to resolve the ambiguity. Inverse 
Compton gamma ray observations provide at least a lower limit to the maximum 
energy of the accelerated electrons, and they provide the number of energetic 
electrons. The latter, in combination with the synchrotron X-ray emission also 
determines the magnetic field strength. 

The signature of inverse compton emission is a high energy continuum. It is 
most important when the radiation field is strong and the plasma density and 
magnetic field are relatively low. 



1.4 Plasma Emission 

Plasma emission is one of the most direct links between microphysical processes in 
a plasma and Earth-based observations. It is most commonly seen in the form of 
type II and type III radio bursts from the Sun. In either case, the radio emission 
is produced by a several step process. A beam of electrons produced in a flare or 
at a shock front penetrates into thermal plasma, giving an unstable bump-on-tail 
velocity distribution. That distribution produces Langmuir waves at the plasma 
frequency as it flattens into a stable distribution. The Langmuir waves can pro- 
duce backscattered Langmuir waves and ion acoustic waves, and subsequent inter- 
actions between the beam-driven waves and these secondary waves produce radio 

1/2 

emission at the plasma frequency of 9np kHz and twice the plasma frequency 



(|Pick fc Vilmeij|2008t ). [Schmidt &i Cairng (j2012l) present an analytical formalism 



for the type II radiatio n from a shoc k. The emission tends to be strongest at nearly 
perpendicular shocks ( Cairnsll201ll) . 



Though plasma emission has been reported from the coronae of a few active 
stars, it is seldom observed from astrophysical sources. Most objects bright enough 
to observe are very dense and optically thick. Type II emission is almost certainly 
produced by shock waves in supernova remants, but radiation at the kHz frequen- 
cies given by the density of the ISM does not reach Earth. 

Since the emission is at the plasma frequency or flrst harmonic, the measured 
frequency directly gives the density in the emitting region. The drift rate of the 
frequency gives the shock speed for a type II burst if the density structure is 
known. However, shock speeds inferred from type II drift rates do not agree well 
with shock speeds measured by coronagraphs i .\lancuxi 20iI7), either because an 
inappropriate density structure is assumed or because different parts of the CME 
shock emit as the CME evolves due, for instance, to selection of a particular angle 
between the field and the shock where emission is efficient. 

The signatures of plasma emission from the solar corona are enormous bright- 
ness temperatures and narrow bands of emission near the plasma fequency. 
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1.5 Dust Emission 

In many cases a plasma is optically thin to radiation from dust, even though the 
emission from individual dust grains is optically thick at some wavelengths. Grains 
in relatively hot plasmas that are heated to temperatures of order 10-100 K emit 
at sub-millimeter and infrared wavelengths. The spectrum is a blackbody modified 
by the opacity of the grain, so it may contain features such as the silicate bump 
at 9.7 microns that can reveal the nature of the grain material. 

Behind the fast shock wave of a supernova remnant, dust is heated to tem- 
peratures around 100 K, mainly by collisions with electrons, even as it is grad- 
ually eroded by sputtering due to collisions with ions. Infrared emission by dust 
can b e the main radiativ e energy loss from shock waves faster than about 300 
km/s (jArendt et al.|[l993 ). The spectrum and the intensity falloff behind the shock 
can be used to infer t he post-shock density and the destruction rate of the dust 
(jWilHams et al.ll2006l[2008ll201ll: ISankrit et al.ll2010h . 

Dust also absorbs and scatters light at optical, UV and X-ray wavelengths. 
The wavelength dependence of the absorption, in particular the 2200 A feature, in 
combination with the IR emission spectrum, is used to infer the size distribution 
and composition of tlie dust (Draine, 2003,) . If the dust column density is fairly 
high, a detectable halo of X-rays appears around a bright X-ray point source 
(|Smith et al.ll2002[) . from which one can derive the grain size distribution and the 
location of the grains along the line of sight. 

The signature of dust emission is a blackbody-like spectrum at IR or sub- 
millimeter wavelengths, sometimes with discrete features due to silicates, poly- 
cyclic aromatic hydrocarbons (PAHS) or other features. The brightness is pro- 
portional to the amount of dust and grain size distribution. It is also sensitive to 
temperature, which in turn is sensitive to either the radiation field that heats the 
grains or the density and temperature of the i gas in which they are immersed. 



1.6 lon-lon collisions 

While most of the radiation detectable at Earth is produced by electrons, en- 
ergetic collisions between ions produce observable gamma rays. These are most 
clearly seen during solar flares, when energetic ions strike the dense gas of the 
chromosphere to produce broad and narrow nuclear de-excitation lines, positrons 
that subsequently annih ilate to produce 0.511 MeV photons, and neutron capture 
lines ( Vilmer et al .l201lh . Gamma ray spectra from RHESSl and other instruments 



can be used to infer the composition of both the chromosphere and the energetic 
ions, the spectral shape of the accelerated particles and their energy content at 
MeV energies. 

At higher energies, cosmic rays can collide with nuclei in the ambient gas to 
produce pious, which can decay into gamma rays. Though it is often difficult to 
tell whether TeV gamma rays are produced by pion decay or inverse Compton 
interaction between ambient photons and energetic electrons, observations of su- 
pernova remnants with the ground-based arrays H.E.S.S, MAGIC and VERITAS, 
and with the FERMI satellite, offer constraints on the acceleration of hadrons 
in stong shock waves. The nature of the gamma ray emission from many SNRs 
is still under debate, but the gamma rays from some old SNRs interacting with 
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dense clouds can be attributed to pion decay fe.g. JUchivama et al ] (l2012l) FERMI 
observations of W44). 



2 Atomic and Molecular Spectral Line Diagnostics 

Atomic spectral lines can appear in emission or absorption. Emission lines usually 
arise following excitation by electron impact or recomb ination into an exc ited 
level, though the y can also be pro duced by ion impact ( Laming et al.l[l996l) or 



photoabsorption ("Noci et al.lll987l) . Cooling by emission of atomic or molecular 
lines often dominates the energy budget of the plasma, and the intensities of 
the spectral lines provide powerful diagnostics for the physical parameters of the 
plasma. This Section provides an overview of atomic emission line diagnostics, and 
Section [3] provides a rigorous discussion of the line formation process. 

Intensity ratios of lines within a single ion can be used to infer the electron 
temperature and density of the gas. Electron temperature diagnostics generally 
hinge on the Boltzmann factor, exp"^^/'^^"^, where AE is the energy difference 
between the two upper levels (Figure [1] left diagram). Such a ratio works best for 
AE ~ kgT, so that optical line ratios are effective for T around 10^ K, where 
AE ~ kgT ~ 1 eV. UV line ratios are effective around lO"' K and X-ray line ratios 
above 10^ K. Often the desirable spectral lines lie at much different wavelengths, 
so that it is hard to obtain a r atio with a single i nstrument, but the technique has 
been applied to solar spectra ( David et al.lll99"8|) . 



The density can be inferred from ratios involving a metastable level. The 
population of that level will be small at low densities. It approaches a constant 
value given by the statistical weight and Boltzmann factor above a critical density 
ncrit = ^2i/ij2i, where A21 is the Einstein A value and 521 is the de-excitation rate 
coefficient. The ratio of a line which involves the metastable level to a line which 
does not will be sensitive to density (Figure [1] center diagram) . Because the A 
values increase rapidly with transition energy and q values decline, Ji^j-it increases 
rapidly from values around 10^ to 10^ cm~^ for optical forbidden lines to 10^ to 
10^° cm~3 for UV hues and lO" to 10^^ cm^^ for X-ray lines. 

A less commonly used density diagnostic takes advantage of the fact that some 
lines formed in the solar corona include both collisionally excited and radiatively 
excited components. The ratio of those components is proportional to the density 
and perhaps plasma velocity and line width (Noci et al. 1987). It is interesting 
to note that for an ion X, the ratios indicated in the center and right panels of 
Figure [1] give < nxn^ / {ricrit + rie) > / < nxne/{ncrit +n-e) > and < nxne > / < 
nxW >, respectively. Here W is the dilution factor of the radiation fSection l3.4p . 
Thus different density estimates are differently weighted averages that do not 
necessarily agree. In principle, comparison of differently weighted averages could 
yield unique information about the distributions of electron density and density 
of the diagnostic ion within t he obs erved volume, but that requires very good 



accuracy for both diagnostics (jLee et a l. 2008.') 



Ratios of emission or absorption lines of different elements can also be used 
to derive the relative elemental abundances. In practice that is often tricky be- 
cause in many cases only 1 or 2 ions of each element can be observed, so the 
ionization state of each element must be accurately known. This usually requires 
a model that involves ionization and recombination rates, each having perhaps 
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Fig. 1 Temperature diagnostics are generally based on the ratio of Boltzmann factors, 
exp~^^/'^B^, in the excitation rates of two spectral lines (left diagram). Density diagnos- 
tics can be based on the competition between radiative decay and coUisional de-excitation 
when the Einstein A value is comparable to the density times the coUisional rate coefficient 
(middle panel) or on the relative contributions of coUisional excitation and radiative excitation 
(right panel). 



a 20% uncertainty, and it often involves an assumption of ionization equilibrium 
(Section 14. 2p that may not be justified (Sections 14.31 and 17. These difhculties 
are somewhat mitigated if one can use ions such as He-like and H-like ions that 
dominate the ionization distribution over broad temperature ranges. 

The profiles of optically thin emission or absorption lines provide a direct 
measurement of the velocity distribution of atoms, molecules or ions along the 
line of sight. Therefore, they provide good diagnostics for the ion kinetic tem- 
peratures, turbulence and in principle non-Maxwellian velocity distributions (Sec- 
tion [5]), though there may be ambiguities among the different interpretations. 

Line profiles directly give ion temperatures when bulk motions do not dom- 
inate. In low density regions of the solar corona, the line widths of oxygen ions 
exceed those of hydr ogen, indicating that the kinetic temperature of O is more than 
16 times that of H (jKohl et al.lll997l: ICranmer et"alll2008l: iFrazin et al1l2003[) . 

Collisionless shock waves are another good example of the application of line 
profile diagnostics. Neutral hydrogen that passes through a strong shock does not 
feel the collisionless shock itself, but finds itself immersed in the hot post-shock 
fiow. Diagno s tics b ased on Balmer line profiles from these shocks are discussed in 
iBvkov et al.l ( 20131) . Most observed line profiles can be fit with a Gaussian or a 
sum of Gaussians, so they are consistent with Maxwellian distributions. The broad 
Ha profile of a bright knot produced by a 2000 km/s shock in Tycho's supernova 
remnant is not Maxwellian, suggesting either a power-law tail or a pickup-ion con- 
tribution, though an interpretation as a sum of Maxwellian contributions cannot 
be excluded (Raymond et al. 201(|). 

Line profiles can be directly used to determine the level of turbulent velocity 
fiuctuations if thermal and bulk velocities do not dominate. Comparison of lines 
from elements of different masses can help to resolve the ambiguity between ther- 
mal and turbulent line widths. Line widths have been used to estima te the level 
of turbulence in reconnection current sheets during solar eruptions ( BemporadI 
Another application has been study of turbulence in interstellar gas us- 
ing the Velocity C oordinate Spectrum me thod to combine line profiles and their 
spatial variations ( Chepurnov et al]|201Cll) . These statistical methods, along with 
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methods based on polarization ( Burkhart et al.ll2012h can reveal the turbulence 
spectrum and whether the turbulence is subsonic. 



3 Optically-Thin Emission Lines 

3.1 Line Formation 

High temperature (> 10^ K) and low density (< 10^^ K) astrophysical plasmas 
are optically-thin to visible, EUV and X-ray radiation. Photons at these wave- 
lengths are generally able to propagate through these environments unhindered 
by opacity effects, such as absorption and re-emission, and scattering, and there- 
fore retain a record of the plasma conditions at the site of emission. Most of the 
radiation in the region of 10® K is due to the emission of photons by electron tran- 
sitions in ions, giving rise to spectral lines. The radiated power per unit volume, 
commonly referred to as the emissivity, depends on: (a) the number of ions that are 
present; and (b) the fraction of those ions in the excite d state that corresponds to 
the transition. For a given transition (in the notation of lMason &: Monsignori Fossil 
11994) : 

P (A,- ,) = Nj A,^,AEj^, [erg cm-^ g-i], (i) 

where Nj [cm "^J is the number density of ions of charge +m in excited 

state j, Aj i [s~^] is the Einstein coefficient and AEj i is the energy of the emit- 
ted photon. The quantity Nj (X^"*) can be rewritten as a series of ratios that 
can be measured observationally or experimentally, or calculated theoretically 
( Mason fc Monsignori Fossilll994) . The total energy flux due to the transition, 
at a distance R from the emitting volume of plasma, can be found by integrating 
the emissivity over the volume and dividing by the surface area of the sphere with 
radius R: 



3.2 The Coronal Model 



A convenient approximation for optically-thin plasmas, such as the solar corona, 
allows a decoupling of the processes that determine the excitation state from those 
that determine the charge state. This can be justified by noting that changes 
in the energy level populations of the emitting ions occur far more frequently 
than changes in the charge state. The processes that determine the excitation 
state are discussed in this Section and those that determine the charge state are 
discussed in Section |31 In optically-thin plasmas energy levels become populated 
by electron coUisional excitation from the ground-state (g) of each ion, and they 
become depopulated by spontaneous radiative decay. It is assumed that timescales 
of photon absorption and electron coUisional de-excitation are far longer. This is 
called the coronal model approximation and in statistical equilibrium the number 
of coUisional transitions from the ground-state g to the excited state j must be 
equal to the number of spontaneous radiative decays back to the ground-state. 
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Ng NeCl^ = Nj A,^g [cm-3 g-i]. (3) 



-,3 



ground-state and level j. If collisions are relatively infrequent then Aj g » NeCg j 
and it follows that Ng >> A^,- (X+"^). There are many more ions in the 

ground-state than in excited states. In a typical transition at EUV wavelengths 
Aj^g = 10^° [s-i] and NeC^j = 1 [s~^] and so for every collisional excitation 
there is an almost immediate radiative decay to satisfy the requirements of sta- 
tistical equilibrium. We note that Equation [3] pertains to 2-level atoms, but ra- 
diative cascades from higher levels, following excitation or recombination, may 
domin ate under particul ar circumstances, such as transitions from Fe XVII 3s 
levels (|Beiersdorfeij|2004 . 



The statistical equilibrium relationship given in Equation [3] and the fact that 

j^(^x+'") ~ ^ leads to an expression for the emissivity in terms of the collisional 
excitation rate: 

N(X+"') n(X) N(H) 
The spectral line intensity is proportional to as expected. 



3.3 Collisional Processes 

The rate at which collisional transitions occur depends on the interaction cross- 
section presented to incident particles by the target and on the flux of incident 
particles. The flux of incident particles can be written: 

F = nvf{E)dE [particles cm"^ s~^], (5) 

where n is the number density of particles, v is the incident particle velocity, E the 
kinetic energy of the incident particles and f{E) the particle distribution function. 
Since particle-particle interactions are mostly via collisions then it is common to 
assume that the distribution function is a collisionally relaxed Maxwellian of the 
form: 

m = ' exp [particles erg-]. (6) 

The electron collisional excitation rate coefficient is found by integrating the elec- 
tron distribution function over the interaction cross-section. 



Clj= / Q^,Jvf{E)dE [cm^s-i]. (7) 

JAE 

AE is the energy difference between level i and j, and this is the lower limit to the 
integral because an incident particle must have at least this much energy in order 
to excite the transition. Qi j [cm^] is the interaction cross-section. In simplified 
form 
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8.63 X 10 



,-6 



Tj J exp 



( 



AE 



) 



(8) 



where ojj is the statistical weight of level i, which is the number of different 
spin and angular momentum states that have energy Ei (the number of degener- 
ate states in energy Ej], arid T^.j (T) is the thermally averaged collision strength 
( Mason fc Monsignori Fossilll994l) . tj^ = 2n^ for hydrogen (where Uq is the princi- 
ple quantum number). 

3.4 Radiative Processes 

Spontaneous radiative decay of electrons from excited states is the dominant de- 
population mechanism in optically-thin plasmas. The generalised radiative decay 
coefficient is: 



The first term of Equation[9]takes account of spontaneous emission. The second 
term accounts for the stimulated component of the emission in the presence of a 
background continuum radiation field, described by a Planck function. In the case 
of radiative decay in the solar corona the background radiation field would have 
a temperature of 5800 K, characteristic of the photosphere. is a dilution factor 
that describes the decay of the radiation field with radial distance, where: 



In the case of the Sun, Rq would be the solar radius and r the distance from the 
centre of the Sun to the height in the atmosphere at which W must be calculated. 
As r — >• CO (e.g. sufficiently far above the surface that r >> Rq) then W — >■ and 
the stimulated component of the emission can be neglected so that Rj i = Aj i. 
The stimulated component of the solar radiative fiux is also negligible at far UV 
and shorter wavelengths; however, photoexcitation of UV lines such as the Lyman 
series and O VI is very important beyond about 1.3 solar radii. 



4 The Charge State of a Plasma 

4.1 Ionization and Recombination 

The charge state of the ions in a plasma is governed by the rate at which elec- 
trons are freed from their bound states and the rate at which free electrons are 
captured into bound states. Bound-free transitions are called ionization and free- 
bound transitions are called recombination. CoUisional excitation and radiative de- 
cay occur on timescales far shorter than ionization and recombination timescales, 
and so these processes can be de-coupled from the excitation and decay processes. 
Ionization (recombination) can then be considered to take place from (to) the 






(10) 
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ground-state of the ion, though it is worth noting that at transition region densi- 
ties (e.g. n ^ 10^" c m~^) ionization and recomb ination from metastable levels can 
become important ( Vernazza fc Ravmond|[l979f ). In optically-thin plasma, such as 



solar and stellar coronae, the important ionization processes are: collisional ioniza- 
tion; and excitation-autoionization. The important recombination processes are: 
radiative recombination; and dielectronic recombination. 

Collisional ionization: as in the case of collisional excitation, the dominant 
process of ionization is by electron collisions (photo-ionization is negligible at the 
energies of interest). Where collisional excitation is generally due to electrons in 
the bulk of the distribution (e.g. a Maxwellian) , ionization arises from electrons in 
the high-energy tail. Since the number density of electrons in the tail is relatively 
low then collisional ionization is relatively infrequent compared with collisional 
excitation. The process of collisional ionization can be w ritten (again employing 
the standard notation o f[M ason 8z Monsignori Fossil 1 1994) : 

X+™-Fe-^X+'"+i-F2e-. (11) 

The ion in the state i loses an electron and a new ion is created in the state i' . 
The incident electron must have sufficient energy to free the bound electron and 
retain enough to remain unbound. 

Radiative recombination: similarly to radiative decay, an important recom- 
bination process is the capture of an energetic free electron into a lower energy, 
bound state, leading to the emission of a photon. The radiative recombination 
process can be written: 



+ e" ^ + AE. 



(12) 



Dielectronic recombi nation: t he dominant recombination mechanism at high 
temperatures, as shown bv lBurges3 ( 1964) . The dielectronic recombination process 
can be written: 



+ e- ^ {X+r) ** ^ X+™ + AE. (13) 

Equation 1131 shows that an ion with m + 1 missing electrons may capture a free 
electron into a particular outer energy level while simultaneously exciting an inner 
electron to a higher energy level instead of emitting a photon. The ( )** notation 
indicates a doubly excited state. The excited inner electron may then decay to 
its original level (or another, if low- lying fine structure states are available) , with 
the emission of a photon, leaving the ion in a singly excited state because the 
captured electron remains in an outer energy level. At this point the recombination 
is complete. Dielectronic recombination is the dominant recombination mechanism 
for most ions at high temperatures, especially those with Anq = transitions 
from the ground state. Dielectronic recombination can also be somewhat density 
dependent, because the emission of a photon often leaves the recombined ion in a 
highly excited state that can be ionized before it decays to the ground state. 

Excitation-autoionization: if the two excited electrons in the second stage of 
Eauation 1131 together have more energy than is needed to remove a single electron 
from the ground state, then the ion is energetically able to autoionize. This means 
that it can decay to the ground state with the ejection of one of the excited 
electrons: 
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(X+™)"^X+"+e-. (14) 

Note that the process described by Equation [Tl] is the inverse process to the first 
stage of dielectronic recombination in Equation 1131 The doubly excited ion has 
two choices: (1) emit a photon; or (2) autoionize (if the total energy of the excited 
electrons exceeds the threshold for ionization). 

Charge transfer: charge transfer between ionized species and neutral hydrogen 
is not usually important in the solar corona, but it can modify the ionization state 
in astrophysical plasmas, especially cool plasma photoionized by a hard radiation 
field. The cross-section for resonant charge transfer is very large and this sometimes 
makes up for a low neutral fraction. 



4.2 The Charge State in Equilibrium 

Ionization and recombination rate coefficients depend strongly on temperature 
and, to a somewhat lesser extent, on density. At higher temperatures the free 
electrons have a greater average kinetic energy and so are able to coUisionally 
release even the strongly bound, inner electrons of the target ions. At lower tem- 
peratures the free electrons are less energetic and can be captured even into the 
low ionization energy, outer bound states of the ions. It is useful to observe that 
ions are typically found at a temperature such that the ionization potential is 
« SfcgT in equilibrium. A full set of ionization and rec ombinatio n rate coefficients 
(e.glArn aud & Rothenfiug 1985; .Arnaud fc Ravmondl fl992: Maz zotta et al.lll998l: 
iBrvans et al.ll2009HDerell2007D for a given element allows the distribution among 
the charge states for the ions of that element to be calculated as a function of 
temperature. We note here that published rate coefficients tend to be calculated 
assuming that the free electrons have relaxed into a Maxwellian electron distribu- 
tion. We consider the consequences of the breakdown of this assumption and the 
calculation and the consequences of departures from an underlying Maxwellian in 
Sections [5l 

One may ask what proportion of helium atoms are neutral, singly ionized and 
doubly ionized at a particular temperature. This is the charge or ionization state 
of the element. At at a temperature of 10^ K helium is fully ionized and so the 
population fractions are: He I (neutral) = 0.0; He II (singly-ionized) = 0.0; and 
He III (doubly- ionized) = 1.0. At 10^ K (adopting the ionization rates of Dere 2009 
and the recombination rates of Mazzotta et al. 1997) the ionization state of helium 
is: He I = 0.0; He II = 0.131; and He III = 0.869. 13% of helium is singly ionized 
and 87% of helium is fully ionized at 10^ K. The population fractions for all the 
ions of a particular element must sum to 1.0 in order to conserve the particle 
number. 

The population fraction for each ion peaks at the temperature at which the 
ionization and recombination rates are equal. More ionizations would act to deplete 
the ion population in favour of a higher charge state, and more recombinations 
would deplete the population in favour of a lower charge state. The ionization states 
for helium given above are only reached when the ionization state is in equilibrium 
with the electron temperature of the plasma. Strictly speaking, as t — > oo at 
T = 10^ K then He I 0.0, He II 0.131 and He HI 0.869. The reason for this 
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is that coUisional processes are not instantaneous. It takes a certain period of time 
for ionization and recombination events to arrange the ions into the charge states 
that correspond to the current electron temperature. As long as the ionization 
and recombination timescales are much shorter than the timescale on which the 
temperature changes then the ionization state can be considered in equilibrium 
with the temperature, and therefore depends only on the temperature. The break- 
down of this condition will be discussed in Section 14.31 

One consequence of de-coupling ionization and recombination from the pro- 
cesses of excitation and radiative decay is that one may assume ionization (recom- 
bination) occurs from (to) the ground state of the ion, and so the rate of change 
of the population fraction of a particular ion i of element X can be written: 

^=n {h-iX,^i + R^X,+i - hX, - i?,_iX,) . (15) 

In the notation of Equation 1151 element X might be helium and then X,_g would 
be neutral helium (He I), and so forth, n [cm^'^] is the electron number density, 
and li and Ri are the temperature dependent total ionization and recombination 
rate coefficients, respectively, with units [cm"^ s^^]. In equilibrium ^ = so that: 

+ = liXi + Ri_iXi. (16) 

The LHS of Equation[T6] comprises the processes that lead to the creation of ion Xi 
(ionization from lower charge states and recombination from higher charge states). 
The RHS comprises the processes that lead to the destruction of Xi (ionization to 
higher charge states and recombination to lower charge states). In equilibrium the 
principle of detailed balance implies that the rate of ionization to Xj is equal to the 
rate of recombination from X^ , and the rate of ionization from X^ is equal to the 
rate of recombination to Xi. This can be expressed in the form of two de-coupled 
equations: 

li-l^i-l = Ri-l^i\ (17) 
Ri^i+l = h^i- (18) 

The ionization state can then be fully specified subject to the final constraint: 

rf=oX, = 1.0, (19) 

where Z is the atomic number of the element X. Making use of Equations 1 1 71 and 
1181 it can be seen that: 

X^-l = ^X, and = ^X,. (20) 

Given a set of ionization and recombination rate coefficients the ionization 
state can be calculated by choosing a suitable value for Xi. The most abundant 
ion i of element X is the one for which Ii{T) ^ Ri{T) at the temperature of 
interest. The population fraction of this ion can then be assigned some arbitrary 
quantity Xi = X'i usually chosen to avoid computational overflow errors since the 
population fractions can vary over many orders of magnitude (this is not so much 
of an issue in the case of double-precision arithmetic) . It is then straightforward to 
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Table 1 Population fractions, and ionization and recombination rate coefficients for the se- 
ries of ions Fe XIV, XV, XVI. These data ar e based on the ionizati o n / re c ombi nation rate 
coefficients from / to a given ion provided bv lMazzotta et"al ] lll998D : lDerj ll2007l) . The rate 
coefficients are in units of [cm^ 

Charge state Population fraction ionization rate Recombination rate 

IQS K IQS K 2.5 X 10*^ K lO" K 2.5 x 10^ K 

Fe XIV 4.60 x 10"" 4.13 x lO^i^ 1.13 x lO^i" 1.35 x IQ-^" 5.16 x 10 

Fe XV 1.41 x 10-5 1.04 x IQ-^^ 6.07 x 10"" 1.04 x lO'^'' 5.09 x 10 

Fe XVI 1.40 X lO"'^ 3.78 x lO"^^ 3.60 x 10"" 9.66 x lO^^^ 2.39 x 10 



calculate [Xl_-^^, X[_2, Xq] and [XIj^-^_, X[^2^ •••i ^'z\ recursively from Eauationl20l 
and find the true population fractions by normalising the values of X[ to 1.0: 

Xi = (21) 



4.3 Non-Equilibrium Charge States 

In circumstances where the electron temperature has been held fixed for a long time 
or the temperature is changing slowly, then the ionization state of the plasma is in 
equilibrium and depends on the temperature only. A slowly changing temperature 
in the present context means that it changes more slowly than the timescales 
on which the processes that change the ionization state of the plasma operate. 
If the temperature change is sufficiently slow then collisions have ample time to 
arrange the charge states of the element such that they are in equilibrium with 
the temperature. 

Consider now a plasma that is heated by some mechanism from 10^ K to 10^ K 
in just 1 second, but it takes several minutes for collisions to change the ionization 
state. In this scenario a plasma of electron temperature 10^ K is created with an 
ionization state (and consequently an emission spectrum) that is characteristic of 
a 10® K plasma in equilibrium. The time derivative in Equation [15] cannot now be 
neglected (however, the bulk velocity will be neglected from the total derivative in 
the following treatment) and a non-equilibrium ionization state arises. In order to 
determine whether non-equilibrium ionization is important in a particular physical 
scenario of interest, the equilibration timescale of the ionization state at the new 
temperature can be estimated from Equation [151 K it is significantly greater than 
the timescale of the temperature change itself then non-equilibrium ionization will 
be important. 

Taking a somewhat less extreme example, suppose that a plasma is heated 
from 10® K to 2.5 x 10® K essentially instantaneously. The equilibrium population 
of Fe XV reaches its maximum at 2.5 x 10® K and so how long does it take to 
equilibrate in this scenario? Based on the data provided in Table [1] we can write: 



Rate of loss of Fe X = n \~Ii{T + AT)X^{T) - Ri-x{T + AT)X^{T)\ 

= nx 1.41 x lO^'^ x (-6.07 x 10"" - 5.16 x 10"") 

= x 1.58 x 10"^^ [s"^]; (22) 
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Rate of gain of Fe X = n[h^i(T + AT)Xi^i{T) + R^{T + AT)Xi+i{T)] 

= n X (1.13 X 10"^° X 4.60 x 10"^ + 5.09 x 10"" x 1.40 x 10 
= nx 5.20x10"" [s"i]; 



Net rate of change of Fe X = n (-1.58 x 10"^^ + 5.20 x 10"") 

= n X 5.04 X 10~" [s-i]. (24) 

For an electron density characteristic of the solar corona n = 10® cm"'^ then the 
equilibration timescale is given by: 

^ ^ 20,000 [s]. (25) 



109 X 5.04 X 10 



If the plasma temperature is changed effectively instantaneously from 10® K to 
2.5 X 10® K and then held constant at the new temperature, then the popula- 
tion of Fe XV will approach equilibrium on an e-folding timescale of 20,000 s. 
In consequence, heating on timescale much shorter than 20,000 s will give rise 
to a non-equilibrium ionization state; for example, the population of Fe XV is 
guaranteed to be out of equilibrium if heating in the solar corona is impulsive (of 
duration shorter than the characteristic cooling timescale) . Heating on timescales 
significantly longer than 20,000 s allows the ionization state to evolve in equilib- 
rium with the electron temperature. We note that even coronal densities of order 
10^^ to 10^^ cm"^ may not be sufficient to maintain the ion population close 
to equilibrium during particularly explosive heating such as occurs during solar 
fiares. The estimate of the timescale provided by Equation 1251 should be 
regarded as an absolute upper-limit. The intermediate population fractions 
of Fe XV and its neighbouring charge states as the system equilibrates are 
not accounted for in the approximation. In essence, the rate of change of 
the population fraction is proportional to the magnitude of the population 
itsel f and it will therefore equilibrate more rapidly as it grows. Figure [5] is 
from lSmith fc HughesI ( 2010l) and shows the characteristic equi libration e- f oldine 



trom lsmitn &: Mughesl l|zUlUl) and siiows trie ciiaracteristic equi libration e- t oldmg 
time-scales for a number of astrophysically abundant elements. iBradshawl (j2009l) 
presents a freely available numerical code that solves the time-dependent ioniza- 
tion equations for all elements up to nickel (Z=28), given any tabulated electron 
temperature and density evolution as a function of time (the time-steps need not 
be uniform). We describe specific examples of scenarios in which non-equilibrium 
ionization might arise in Section [3 



5 Non-Maxwellian Electron Distributions 

Temperature changes on timescales much shorter than those on which ionization 
and recombination can change the charge state of the plasma are not the only way 
in which the ionization state may be different than expected for a given temper- 
ature. If the electron distribution function is driven away from Maxwellian with 
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Temperature (K) 



Fig. 2 'Smith & Hughes' ('20101): the left axis measures the density- weighted timescale [cm ^ s] 
for several abundant elements to achieve one e-folding toward ionization equilibrium in a 
constant temperature plasma; the right axis measures the density-weighted timescale for the 
elements to reach within 10% of their equilibrium population. 



the addition of a significant population to tfie higti-energy tail of the distribu- 
tion, then ions of greater charge can be created at some temperature that is lower 
than the temperature at which they arise in equilibrium. Non-Maxwellian distri- 
butions can arise in several ways; for example, in an astrophysical context they 
may be expected to occur in circumstances where a region of very hot plasma is 
separated from a much cooler region by a steep temperature gradient, as is the 
situation in the solar atmosphere. CoUisionless electrons may then stream from 
the hot, less-dense plasma down the temperature gradient into the cooler, denser 
plasma driving the tails of the electron distributions in these regions away from 
Maxwellian. Non-Maxwellian distributions can also be induced when a beam of 
particles is accelerated by some mechanism, such as magnetic reconnection, and 
interacts with the background plasma. A non-Maxwellian distribution can be cre- 
ated in laboratory plasma by laser-heating. 

Departures of the electron population from Maxwellian distributions have im- 
plications for several properties of the plasma, among them the excitation and 
ionization states (and, consequently the spectral emission) of its component ions 
and the transport of energy by thermal conduction. It is therefore highly desirable 
to take account of these effects in numerical modeling studies, but this is an ex- 
tremely difficult task to achieve self-consistently. There are two general approaches: 
(1) carry out particle-in-cell type calculations where the distribution functions can 
be obtained directly; (2) carry out calculations based on the fluid equations de- 
rived by taking successive moments of the underlying distribution function. The 
first approach is discussed elsewhere in this volume (REFERENCES TO AR- 
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TICLES IN THIS VOLUME TO BE ADDED). The difficulty of the second 
approach is that solutions to the fluid equations can quickly become inconsistent 
with the assumptions on which their derivation is based. For example, in the case 
of steep temperature gradients the mean-free-path of even thermal electrons can 
approach (and exceed!) the temperature scale length and then the plasma cannot 
be considered collisional on the characteristic scale length of the fluid. However, 
the validity of the fluid equations is contingent on the coUisionality of the sys- 
tem on the relevant spatial scale. One advantage of the second approach over the 
first is that plasma systems can be modelled across a much larger range of spa- 
tial scales. In particle-type codes one is generally confined to studying phenomena 
on a particular scale, such as the width of a conduction front or a shock, or the 
scale of the diffusion region in reconnection (REFERENCES TO ARTICLES 
IN THIS VOLUME TO BE ADDED). In fluid codes, the solution can range 
across many spatial scales from a few meters to hundreds of thousands of kilo- 
meters (e.g. in the case of the Sun's atmosphere). The challenge is to develop a 
method by which physical phenomena on particle scales can be self-consistently 
included, when needed, in a code that operates predominantly on fluid scales. This 
requires the distribution function to be calculated from a suitable kinetic equation 
in tandem with the time-advancement of the system of fluid equations in order 
that corrections can be made to the fluid variables. 



5.1 Kinetic Equations 

The key to calculating the distribution function in a collisional or weakly-coUisional 
plasma is the manner in which the collision term of the Boltzmann equation is 
treated. One of two approaches is usually adopted. The flrst is to handle collisions 
via a phenomenological term based on the expectation that the plasma particles 
will ultimately relax into a collisional / Maxwellian distribution on some timescale 
th at depends on the degree of coUisionality of the plasma. This was first suggested 
bv iBhatnagad ( 1954|) and kinetic equations of this form are commonly referred to 



as BGK equations (based on the initials of the authors of that paper). 



(f) 



, „. . = l^SS {Fs - fs) + I^SS' {F,,, - fs) , (26) 

C£ / collisions 

where s,s' denote the particle species (e.g. electrons and protons), Vssjfss' are 
the species and inter-species collision frequencies, and F denotes a Maxwellian 
distribution determine d by the loca l properties of the plasma (e.g. temperature, 
density and bulk flow) . iMorsd ( 19631) studied the energy and momentum exchange 
between non-equipartiti on gas es in the cases of Maxwell, Coulomb and hard sphere 
interactions, and , Morsel (|l964j) showed how to choose free parameters for the cross- 
collision terms in BGK-type models to conserve density, momentum and energy. 
This work was limited by the underlying assumption that Ue/rei = ni/ri^, which 
for a fully-ionized hydrogen plasma {jie = n.i) implies that electrons and ions are 
equally affected by their mutual collisions (rgj = Tje) when in reality they relax 
on a tim escale longer by the square root of the mass ratio (r^e = ^JrOifrneTi^i). 
iGreend ([1973) then developed a simple improvement for BGK-type models of 
electron-ion collisions to produce the correct relation between the time scales of 
ion-electron momentum exchange and ion thermalisation. This work showed how 
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to obtain the correct timescale ordering and how to choose the correct parameters 
for the Maxwelhans in the cross-colhsion terms to conserve density, momentum 
and energy. 

The second approach to handling collisions is to assume that changes in the 
velocities of charged particles are due to the cumula tive effect of long-range en- 
counters via inverse square forces (e.g. Landaul 19361) . The collision integral can 
then be written (e.g. iLiepoievic &l Burgesj 1990bl) : 



=-S,— (fs(AvlS\ +\E,,^^(fslAvlAvi)'\, (27) 
V 9t /collisions 'dvlV \ I) 2 dv\dvl\ \ 1) ^ ' 

where 



Av\ 



fs' (^ 



' {g, n)gAvld^nd^Vgi 



(28) 



and 



(^AvlAvi^^ = E,, J f,, {V,,) a,,,{g,n)gAvlAvid^ad\',,. (29) 

The ter ms of Equations[27]to[29]are described in detail in Section 2(a) of lLiepoievic &: BurgessI 
(|l990bj) . Equati ons of the form of [ 271 are commonly referred to as Fokker-Planck 
(FP) equations. ICohen et al ] (|l95(]|r adopted a method of approximating the dis- 
tribution function by representing it as a Maxwellian plus a small perturbation to 
calculate the electrical conductivity of a gas. Their approach is valid in the pres- 
ence of weak spatial gradients and weak electromagnetic fields. The distribution 
function then takes the form /s = /o + fi where /p = Fs and 



/i = FsD{vs)fi, 



(30) 



where fi is the cosine of the pitch angle. D is a function of the particle speed 
found by substituting /s = /o + fi for fs in the Boltzmann equation, linearis- 
ing the co llisional integral in /i, and solving the integro-differential equation. 
ICohen et~al . ( 1950) n eglect e d elec tron-electron interactions from the collision op- 
erator but Spitzer fc HarmI ( 19531) . in what is now considered the 'classical' treat- 
ment, followed the same approach and included electron-electron interactions in 
their collision operator. They also extended the solutions to completely ionized 
gases with different mean nuclear charges and calculate the electrical and thermal 
conductivities of the gas. The solution to the integro-differential equation in the 
classical treatment has the form (for electrons) 



fe=Fe[l~X 



eE 



ZDe ( 

A \kBTe 



1 



dPe 

ds 



^ ZDt 1 dTe 

' B T^lh 



(31) 



The quantites ZD^/A and ZDq^/B are tabulated in [spit zer fc HarmI (|l953l) as 
functions of the electron speed normalised to the thermal speed and Aq is the mean- 
free-path of thermal electrons. Above a particular speed Uj^j-j^ the calculated value 
of /i becomes comparable to /o and t he assumptions made to derive Equation [31] 
are invalidated. The approximation of ISpitzer fc Har"ml ( 19531) is only valid in the 
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low-velocity regime v < v^,^.^^ and the upper limit of the regime depends on the 
str ength of the electri c field E, and the temperature and pressure gradients. 

iRosenbluthI ( 19571) derived the Fokker-Planck equation for arbitrary distribu- 
tion functions in the case where two-body interactions are governed by a force 
that obeys the inverse square law. The coefficients Av and AvAv in the Fokker- 
Planck operator were written in terms of two fundamental integrals / potentials 
that depend on the distribution function of the background particles (including 
those of the same species) . Expanding the distribution function as a set of Legen- 
dre functions of the pitch angle, the Fokker-Planck equation is cast into the form 
of an infinite set of one-dimensional, coupled non-linear integro-differential equa- 
tions. Approximating the distribution function by a finite series, the Fokker-Planck 
equations can be solved num erically. Keeping one te rm of the series corresponds to 
the approxima te solution of |Chandrasekhar'V l943) an d keeping two terms yields 



the solution of I Cohen et al.l (,195 0). RosenblutbTi 19571 ) showed that 



( 



dfs 
dt 



-r- 



' collisions 
with the Rosenbluth potentials 



dvl 



dvl 



1 d' 



2 dvldvi \ dvldvi 



f- 



rG 



(32) 



H(vs) = Es> 



ms + m. 



d^y's'fs' (Vg') \vs 



where F : 



4TT{Zsef{Zsef\iiAs 



(33) 



iLiepoievic fc Burges^ (|l990bl) presented a step-by-step description of a method 
for calculating the distribution function in the presence of strong gradients. In 
their method th e low-velocity part of the distribution is given by solution of 
zer fc HarmI (jT953). The high- velocity tail of the distribution function is given 
as a solution to the high-velocity form of the Fokker-Planck equation which is 
derived from Equations 1321 and 1331 bv neglecting the interaction between the high- 
velocity particles themselves and considering only their interaction with the low- 
velocity (near Maxwellian) part of the distribution function. One may then derive 
a linearised form of the Boltzmann equation with the Fokker-Planck collision op- 
erator that applies to high-velocity particles. For electrons: 



(t) 



collisions 



V, 



1 d 

^ dVe 
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Vf^u{Ve 



kpTe dfe 
dVe 



me 



+ Veft 



/i 2\ dfe 



The full derivation of Equation [3l] is given by IlI 



lepoievic . 



, (34) 

Burgess (l990b) on 

pages 73 to 88 of their article. They also describe in detail a numerical treatment for 
its solution following non-dimensionalisation and transformation into a form more 
convenient for numerical work. The solutions in the low- velocity and high- velocity 
regime are combined, subject to suitable matching conditions (e.g. smoothness), 
in a region of the velocity space where both methods are approximately valid; two 
thermal speeds was found to be the optimal value. 
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Given the significant complexities that are involved in working with the Fokker- 
Planck equation it is tempting to revert to B GK-type approximatio ns of the colli- 
sion operator. However, one must be careful. [Livi fc MarschI ( 19861) compared the 
coUisional relaxation of a double-beam and a bi-Maxwellian distribution function 
for a Fokker-Planck and a BGK collision operator. They found that moments of the 
distribution function up to and including temperature (the 2nd moment) were in 
good agreement between the two schemes when the frictional energy-loss rate was 
used as the effective collision frequency in the BGK operator, but that the heat flux 
(the 3rd moment) exhibited differences due to its sensitivity to the shape / skew 
of the distribution function, which enters the Fokker-Planck operator via the sec- 
ond d erivative of the distribution function w.r.t. velocity. Liepoievic & MacN eicd 
( 1988t ) calculated contributions to the heat flux in a solar flare atmosphere from the 
tail of the distribution functi on using the high- velocity form of the Fokker-Planck 
equation (jLiepoievic fc Burg ess 1990b) and compared the results with a BGK- 
type calculation. They found that the BGK technique can estimate contributions 
from the high-energy tail to the heat flux to within order of magnitude. 

As computers became more powerf ul, deta i led n umerical treatments of the 
Fokker-Planck equation became feasible. Shoub (19831) provides a detailed discus- 
sion of the break-down of the ,Spitzer fc Harm (1953i) calculation of the electron 
distribution function and describes an approach to deriving and then solving nu- 
merically the high-velocity form of the Landau-Fokker-Planck equation. Implica- 
tions of the break-down of the local Maxwellian approximation are discussed for: 
energy balance in the upper chromosphere and low TR; the He resonance line spec- 
trum; the Schmahl-Orrall observation of continuum abs orption by neu tral H, and 
the origin of the 20,000 K temperature plateau. However. lShoubl ( 19831) was unable 
to say anything quantitative about the heat flux since the kinetic equation was only 
solved to six thermal speeds. Had IShoubl (|l983l) a pplied the same transformation 
following non-dimensionalisation as employed by iLiepoievic fc Burgesd (|l990b|) 
then it would hav e been p ossible to significantly extend the calculation in velocity 



space. 



Liepoievi 3 (|l990a|) used the approach described in ILiepoievic fc Burgessl 



()l990bl) to show that dist ribution functions are near Maxwellian in the commonly 
used FAL ( Fontenlalll99 3') model s of the ph otosphere to mid-TR and so the models 
are valid in their given form. MacNeice et a l. (1991) applied the same approach to 
the transition region of a flaring loop and found a substantial enhancement in the 
tail populations throughout that region of the atmosphere. We discuss the results 
of some attempts to apply these calculations of distribution functions to fluid mod- 
els in order to take account of the consequences of non-Maxwellian distributions 
in the following Section. 



5.2 Heat Flux / Transport 

The fluid equations are derived by taking successive moments of the Boltzmann 
equation when it is written in terms of distribution functions that exhibit only 
small deviations from a fully-relaxed Maxwellian distribution. Since the statistical 
treatment of a particle ensemble in terms of a fluid is valid only in this coUisional 
limit, then only small deviations can be tolerated. In general, departures from 
Maxwellian are treated as a perturbation and the distribution is expanded in terms 
of some parameter that should remain small, such as the ratio of the electron mean- 
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free-path to the temperature scale length (the Knudsen number, Kn), in order to 
derive non-linear terms of the fluid equations such as the heat flux. 

The transport of heat by thermal conduction is the dominant transport process 
in hot but tenuous astrophysical plasmas. It determines the temperature and thus 
the density structure, via the temperature-dependent scale length, in the solar at- 
mosphere (for example) and so it is important to handle it as accur ately as possible. 
The m ost commonly used form for the heat flux is that given by ISpitzer fc HarmI 
(|l953h . valid in the limit of weak gradients and weak electric fields: 



Fc = -kVT, (35) 

where the conductivity k = kqT^^"^ (for a fully ionized hydrogen plasma) and the 
constant is the quantity calculated by Spitzer {kq ^ 10~®). Despite its strong 
non-linearity the form represented by Equation 1351 is convenient to implement in 
fluid-based numerical codes, but cannot guarantee an accurate representation of 
the heat flux if it is used indiscriminately. Experimental and numerical results 
have shown that its range of applicability is actually quite limited. Equation 1351 
indicates that the heat flux can increase indefinitely provided that the tempera- 
ture gradient continues to steepen, but eventually a physical limit must be reached 
when there are no more particles remaining to support the implied heat flux. This 
is the free-streaming limit, essentially the maximum heat flux that the plasma 
can sustain, and may be estimated by assuming that the majority of the parti- 
cles (e.g. electrons) stream down the temperature gradient at the thermal speed 
( Bradshaw fc Cargillll2006h (more sophisticated numerical treatments indicate the 



free-streaming limit is about 1/6 of this value). At the very least, then, a limiter 
should be deployed in any numerical model that uses Equation 1351 in order to 
constrain the heat flux to physically justiflable values. 

There have been a number of efforts to derive systems of fluid equations that 
take account of stronger depart ures from Maxwellian to be implemented in nu- 



merical models. ICampbell ( 1984h found a solution to the Boltzmann equation that 



extends the Chapman-Enskog approximation to large temperature gradients and 
electric flelds, to calculate electron transport in a fully ionized gas. The collision 
term was written in the form of a collisional relaxation with a velocity-dependent 
relaxation time deflned in terms of the scattering length. The distribution function 
was assumed to be separable with the angular dependence represented by a slowly 
varying function. Calculating th e moments of this dist ribution function led to cor- 
rection factors to the classical ( Spitzer fc Harmlll953h transport coefficients as a 



function o f the tempe rature gradient scale- length and an inherently flux- limited 
heat flow. iKillie et al.l (j2004h derived a complete set of fluid equations for fully 
ionized gases that improve the treatment of Coulomb collisions by taking into ac- 
count the shape of the distribution function to better calculate the heat flux and 
the thermal force. They chose an analytical velocity distribution function with a 
Maxwellian core plus a high- velocity correction term proportional to v'^, and ob- 
tained transport equations by inserting their choice of distributio n function into 
the Boltzmann equation with a Fokker-Planck collision operator. IChiuderi et al.l 
(l2on) derived a set of two-fluid equations applicable to weakly collisional plas- 
mas by using a relaxation approach to the collision operator and selecting 'mixed' 
Maxwellian distributions for the two interacting species that conserve momentum 
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and energy. The collisional term in their treatment depends on an 'average' or 
're presentative' collis i onal t imescale that is velocity- independent. 

iGrav fc Kilkennvl ( 1980f ) described the results of experiments in which the ra- 
tio of the electron mean-free-path to the temperature scale-length was found to 
be about 0.5 and enhanced low-frequency turbulence was observed. They used a 
numerical simulation of the experimental set-up to show that ratios of 0.5 implied 
a heat flux limited to less than 5% of the free-streaming limit in the hot part of 
the plasma. They also found Te/T^ in the same region sufficient to excite heat- 
flux driven ion acoustic turbulence, thus explaining the low frequency turbulence 
observed in the experiment. The observed level of turbulence in the experiment 
was enough to account for the predicted low thermal conductivity in the numer- 
ical model, which w as due to e lectro n sc attering from interactions with the ion 
acoustic turbulence. iBeh et al.l (|l98ll) and lMatte fc Virmond (|l982l) studied elec- 
tron heat transport down steep temperature gradients in laser-induced plasmas 
by numerically solving the Fokker-Planck equation. The heat flux was found to 
be substantially smaller than that predicted by the classical theory or the free- 
streaming value when the mean-free-path reac hed a fraction of only one-hundredth 
of the temperature scale length {Kn = 10~^). [Shoubl ( 19831) found significant de- 
viations from Maxwellian in the tail of the distribution for Kn = 10~^, but was 
unabl e to provide a quantitative estimate of the heat flux. lOwocki fc Canfieldl 
( 1986[) used a BGK-type method to calculate the electron distribution function 
in the solar transition region to study the effect of a high-energy tail on the heat 
transport and collisional excitation and ionization rates. For the case studied they 
found that non-classical transport does not significantly alter the excitation or 
ionization state of ions with emission lines that form predominantly in the lower 
transition region (with excitation energies in the range 10 eV, because electrons 
with these low energies thermalise quickly), but the non-classical heat fiux in this 
region does depend sensitively on the temperature gradient in the upper transition 
region. 

In the case of pronounced departures from Maxwellian distributions it is clear 
that correction factors and localised approaches to calculating the distribution 
function, and hence the heat flux, are not sufficient. For example, contributions to 
local quantities from non-local sources may lead to strong departures from local 
Maxwellian distributions. Such kinetic behaviour is inherently incompatible with 
the fiuid approximation in which it is assumed that the properties of the plasma 
can be determined entirely locally (e.g. the heat fiux as a function of temperature 
and the temperature gradient). Since it is generally not feasible to solve a kinetic 
equation (certainly not a time-dependent form) in tandem with the fiuid equations 
to correct for the consequences of kinetic behaviour, then the challenge is to find 
an alternative; e.g. a computationally tractable approach that can be implemented 
in an otherwise fiuid-based treatment, and that permits one to account for purely 
kinetic effects (e.g. non-local infiuences) on quantities such as the heat fiux in 
regions where the Knudsen number grows large. One such approach is to adopt a 
del ocalisation formula. 

iLuciani et al. 1 (|l983l) found a delocalisation formula for the heat fiux, using a 
set of solutions to the Fokker-Planck equation. Delocalisation formulae are based 
on delocalisation kernels that operate on calculations of the heat fiux made using 
Equation 1351 The kernel essentially acts to 'smear' the classical heat fiux out 
over the computational domain in a manner that mimics the spatial profile of the 
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heat flux tha t would be found fro m a full Fokker-Planck calculation. The formula 
presented bv lLuciani et alJ ( 1983h has the form: 



and the delocalisation kernel is 



Fc{s)= / w{s,s')FsH{s')d.s' 



(36) 



d(s,s') 



2A(s') 



exp 



X{s')n{s') 



ds ' 



(37) 



The quantity A is an effective range for the electrons, related to the mean- 
free-path. In the limit of shallow temperature and density gradients the kernel 
w behaves like a 5-function, where J ■w{s,s')ds' = 1 and Equation 1361 reduces to 
Fc = Fgjj. Despite the double integration, the delocalisation formula is straight- 
forward to effici ently implement in a fluid code to replace the heat flux in the form 
of Equation l35llLuciani et aL (Il985 ) found an analytical justification for the delo- 
calisation formula and Bendib et al. ] (|l988l) developed an impr ovement that takes 
the pr esence of an electric field into account. More recently, lAlouani-Bibi et aD 
(|2004h studied non-local electron heat transport using a number of different ap- 
proximations to the Rosenbluth potentials in the Fokker-Planck equation to find 
delocalisation kernels for non-local heat flux formulae to be used in fluid codes. 

A number of authors have compared the different approaches to calculating the 
heat flux and have implemented them in numerical models in order to apply them 
to particular problems in which non-Maxwellian electron distributions are expected 
to arise. [Smith (1986) discussed classical (Spitzer & Harm 1953), locally limited 
(|CampbeIllll98l and non-local ( rLuciani et aLHl983Hl985t) algorithms for the heat 
flux and their application to heat transport in the case of the steep temperature 
gradie nts (thin conduction fron ts) that arise during the impulsive phase of solar 
flares. iKarpen fc DeVor^ (jl987|) investigated how these different formulations for 
the heat flux affect the physical characteristics of the corona, transition region and 
chromosphere in numerical models of solar flares. Both sets of authors found that 
the heat flux in the hot part of the plasma obtained with the non-local treatment 
was smaller than the locally limited and classical values, whereas the heat flux in 
the colder parts of the plasma (e.g. in the transition region and chromospere) was 
signiflcantly enhanced compared with the locally limited and classical values. In 
consequence, both flux limiting and delocalisation play an important role in the 
evolution of the plasma. In the case of flares this leads to a 'bottling up' of energy 
in the corona, allowing it to reach much higher temperatures, and the earlier onset 
of weaker chr omospheric evaporation. 

iLiepoievic fc MacNeicd (jl989) calculated the heat flux in a solar active region 
coronal lo op from distribution functions obtained using the sophisticated model de- 
scribed in lLiepoievic fc Bureesd (Il990bl). and compa red it with the heat flux given 
by the classical treatment of Spitzer fc Harm (1953j, the correction coefficients to 
the classical treatment gi ven by [ Campbell (1984^') an d the heat flux given by the 
delocalisation formula of iLuciani et al. ( 19831 1985|) (Figure [3} . They concluded 
that the classical treatment failed completely in the lower corona, predicting a 
strong heat flux flowing down the temperature gradient when the kinetic equation 
yielded heat flux flowing up the temperature gradient, and the possibility that the 
role of the heat flux could be misinterpreted in the energy balance of the corona. 
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Fig. 3 The heat flux in a s o lar act ive region loop calculated by iLiepoiev ic &: MacNeicj 
||1989|): iLiepoievic fc Burgesd jl990bl1 (solid hne),\.Spitzcr ,_fc Harm (1953) (dotted line), 
Campb(d7fl98^) (dashe d line) , and llLuciani et"al IllQSlllQSSl) (dot-dashed line). Image credit: 
Liepoievic fc MacNeicS ||1989I) . 



iLandi fc Pantellinil ( 200lj) also found that the heat flux can flow up the temper- 
ature gradient in the case of supra-therm al tails characterised by k distributions 
with K < 5. Encouragingly. ILiepoievic fc M acNeice (1989) did find relatively good 
agreement between the deloca lisation formula and the more sophisticated kinetic 
calculation. IWest et al.l (l2008|) implemented the delocalisation formula in the HY- 



DRAD (e.g. iBradshaw fc Masoiil 120031: iBradshaw et al.ll2012l) code to investigate 



the lifetime of hot, nanoflare-heated plasma in the solar corona. The aim of this 
work was to determine whether the bottling up of energy in the corona due to the 
severe heat flux limiting that arises in the limit of large Knudsen numbers pro- 
vided s ufficient time for the ionization state to equilibrate following rapid heating. 
lAlouani-Bibi fc Matte (2002 ) developed a non-local model of electron heat flow 
in laser-heated plasmas, taki ng into account super-Gaussi an deformation of the 
electron distribution function. lAlouani- Bibi fc Matt^ (|2003l) derived an analytical 
description of electron-ion energy exchange by Coulomb collisions in the presence 
of super-Gaussian electron distributions, and found the ratio Ti/Te at which the 
coUisional electron-ion energy exchange cancels increases from 1 in a Maxwellian 
plasma to 1.98 in a super-Gaussian plasma. 



5.3 Excitation, Ionization and Radiation 

The specific nature of the local distribution function can have an important effect 
on the rate of excitation and ionization via collisions. Excitation is generally a 
consequence of interactions between ions and electrons in the bulk of the distribu- 
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tion, but ionization is particularly sensitive to the tail population. The stronger 
heat fluxes at the base of steep temperature gradients found in a number of the 
studies described in Section [5 . 21 imply enhanced tail populations of streaming elec- 
trons, which can feasibly alter the ionization state such that it can no longer be 
considered a strong function of the local temperature (and, to a lesser extent, the 
local density) alone. CoUisional excitation rate coefficients can be calculated by 
substituting a suitable distribution function into Equation [7] and the ionization 
rate can be calculated by inserting the appropriate ionization cross-section (usu- 
ally pertaining to the ground-state) in place of Qij- Investigations of the effect 
of non-Maxwellian distributions on the ionization state have proceeded along two 
general lines; (a) calculate the distribution function by solving some simplified 
form of the Boltzmann equation (e.g. BGK, Fokker-Planck) ; or (b) choose an an- 
alytical form for the distribution function with the properties of a MaxweUian at 
low-velocities / energies, but which permits an enhanced tail population where 
the degree of enhancement can be controlled by a single parameter. A popular 
generalisation of the MaxweUian distribution that fulfils these requirements is the 
K-distribution: 



3/2 



(K-1.5)kBT 

A. = 3,, . (39) 

r(K-0.5)(K-1.5)^'^^ 

The K-distribution has the form of a MaxweUian in the limit k — )■ cx). The most 
probable energy of a particle in the distribution is Ep = {k — l.5)kBT/K and the 
mean energy of the distribution is < _E >= 3kBT/2 (i.e . independe n t of k and 
the same a s the MaxweUian at the same temperature). lYoon et al.l (j2006l) and 
iRhee et al.l (|2006l) demonstrated that K-distributions can be induced by sponta- 



neous scattering (absent in collisional treatments) when electron beams are accel- 
era ted by weakly turbulent p rocesses . 

lOwocki fc Scuddeil (|l982l) used ^-distributions to study the ionization state of 
gases with non-Maxwellian electron distributions, finding changes from the ion- 
ization temperature assuming an underlying MaxweUian distribution of up to a 
factor of 2. The importance of the high- velocity tail to the ionization state de- 
pends on the ratio of the ioni z ation potential to the mean thermal energy of the 
electrons. lOwocki fc Scuddeil (|l982f) also found that the high ionization energy 
required for the O VIII -f-*- O IX transition means that oxygen ionization at solar 
coronal temperatures is more sensitive to the tail of the distribution than ele- 
ments of lower ionization ene r gy (su ch as iron) found within that temperature 
range. iDzifcakova &: Kulinoval (|2003l) calculated excitation, ionization and recom- 



bination rates for ^-distributions for a range of values of the parameter k. They 
found changes in the level populations and the relative ion abundances. A syn- 
thetic spectrum was calculated which showed that some C III, C IV and O IV 
lines are sensitive to the shape of the distribution function and their intensities 
enhanced by a factor 2 — 6 in the presence of strongly non-thermal distributions. 
IPzifcakova (2006) investigated the influence of K-distributions in the solar corona 
on Fe VIII - XV excitation and ionization, and on the line intensities associated 
with those transitions. They concluded that it ought to be possible to diagnose 
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Fig. 4 The X-ray and UV continuum at 10 MK calculated from k and n distributions rep re- 
senting different strengths of departure from Maxwellian. Image credit: Ibudfk et al.l ||2012| ). 



the value of k that would best characterise the electron distribution from ratios 
of Fe IX 171 A, Fe XII 195 A and Fe XV 284 A lines provided that the plasma 
density is known. 



iDzifcakova fc MasonI (|2008ah calculated non-Maxwellian electron excitation 
rates for ions of astrophysical interest. They demonstrated a method for extracting 
collision strengths from the Maxwell-averaged collision strengths (Upsilons) that 
are provided by the CHIANTI atomic database and then integrated these over 
the specific non-Maxwellian distribution in order to calculate the corresponding 
excitation rate. K-distributions, employing a range of values of k, were used to cal- 
culate synthetic spectra for Fe XV and XVI in the 50 — 80 A range for comparison 
with solar observations. For c onsistency in the generation of the synthetic spectra, 
IDzifcakova fc Mason (l2008a ) used t he equilibrium ionization states of Fe derived 
for a range of k by Dzifcakovil ( 20021) . However, they found no conclusive evidence 
for non-Maxwellian distributions in the particular flare dataset that was compared 
with the synthetic spectra. 



iDudfk et al.l (j201lh calculated the bound-bound and free-free radiative losses 
arising from plasmas with non-Maxwellian electron distributio n functions using k- 
and n-distributions (Equations 5 and 6 of iDudi'k et al.ll201lh . It was found that 
changes in the radiative loss function due to non-Maxwellian distributions are 
greater than errors in the atomic data and errors due to missing contributions 
from free-bound continuum. While radiative loss functions for k distributions are 
generally weaker than for Maxwellians, the opposite is true for n-distributions. 
They also found that the contribution from bremsstrahlung changes by only a few 
percent, except in the extreme case of k = 2. Following on from this earlier work, 
iDudik et al.l ( 2012|) calculated the X-ray, UV and radio continuum arising from 
non-Maxwellian distributions using k- and n-distributions (Figure H}. They found 
that at flare temperatures and hard X-ray energies both the bremsstrahlung and 
the free-bound spectra are dependent on the assumed distribution, and concluded 
that the low energy part of k distributions can be determined from observations 
of the continuum. 
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5 The optically-thin radiative loss function. Image credit: [Landi &: Landinil 11199911 . 



6 The Radiative Loss Function 



The power per unit volume of plasma emitted by a single spectral line is given by 
Equation m The total power per unit volume is then the sum total of the power 
emitted by the many thousands of spectral lines which belong to the chemical 
elements that comprise the plasma. In the case of optically-thin astrophysical 
plasmas this quantity can be written in a conveniently compact form: 

Er = NeNnAiTe) [erg cm-^ s"!], (40) 

where N^j is the number density of hydrogen atoms. In a fully ionized hydrogen 
plasma A^'e = Nfj. En is generally referred to as the radiative volumetric loss rate. 
A(Te) is called the optically-thin radiative loss function (or the total emissivity of the 
plasma, as shown in Figure [S]) and it encapsulates a vast amount of atomic data. 
The radiative loss function depends upon the element abundances, the ionization 
state of the elements, and transition energies and probabilities. These must be 
determined for hundreds of ions and in many cases thousands of spectral lines per 
ion, in order that accurate radiative loss functions can be calculated. As atomic 
data is improved and updated then so must the radiative loss function. The most 
convenient way to keep abreast of developments is to use a comprehensive and 
regularly updated atomic database, such as Chianti ( Dere et al.lll997i : lLandi et al.l 
[2012. : https:// www.chiantidatabase.org) , which contains carefully assessed atomic 
data and the appropriate functionality for calculating spectra and A. A number of 
assumptions concerning the properties of the plasma, such as the nature electron 
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distribution function and the time-dependence of the ionization state, are implicit 
in the most commonly used calculations of A for astrophysical applications. The 
limits of these assumptions have been demonstrated in Sections [3] and [S] and so 
the emissivity of individual spectral lines and the radiative loss function must 
be recalculated, whenever these limits are reached, to take proper account of the 
microphysical processes that can arise in astrophysical plasmas. 

Following Equation[3]the emissivity of a single transition between energy levels 
j and g m a. particular ion is given by 

_ N {X+"^) N{X) N{H) ^ergcm3s-M (41) 

N{X) N{H) Ne ^1,0^^3,9 lergcm s J. (41] 

The ion emissivity is then obtained by summing over all of the transitions that 
may occur within the ion: 



Ax^ = ^Aej, g, (42) 

so that Axi is the radiative loss function for the particular ionization state i of 
element X. The radiative loss function for the element can be found by summing 
over the number of charge states: 



Ax = S,Ax, (43) 

and the total radiative loss function is the sum over the number of elements of 
interest: 



A = SxAx- (44) 

The radiative loss function is strongly dependent on the electron temperature Te 
in equilibrium, but it is clear from Equation Hi] how important the ionization 
state (the first factor on the right-hand side) and the coUisional excitation rate 
Cgj are to accurately calculating it. When the ionization state exhibits strong 
departures from equilibrium then the temperature dependence of the radiative 
loss function can be lifted, the rate of energy loss by radiation, and the intensities 
of individual spectral lines, may not be characteristic of the actual electron tem- 
perature. Furthermore, the ionization state (via the ionization and recombination 
rate coefficients) and the collisional excitation rate depend on the underlying elec- 
tron distribution which is generally assumed to be Maxwellian, but circumstances 
can easily arise astrophysical plasmas when this assumption is certainly not valid 
(Section O. 



7 Signatures and Diagnostics of Non-Equilibrium Processes 

Analytical analyses can identify the conditions under which non-equilibrium pro- 
cesses become important to understanding the properties and behaviour of as- 
trophysical plasmas, and numerical models demonstrate that such conditions are 
commonplace in the optically-thin astrophysical plasma systems that are the focus 
of a great deal of current research interest. In this Section we consider potential 
signatures of non-equilibrium processes and the evidence for their manifestation 
in observational datasets. 
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7.1 Non-equilibrium Ionization 



iGriemI ( 1964|) discussed the potential for departures from equilibrium of the ion- 



ization state in cases where the dynamical evolution of the plasma occurs on 
timescales that are shorter than those of ionization and recombination. He cited 
the particular example of transporting the ion population across a st rong temper- 
ature gra dient, as might be the case in the solar transition region. Ijoslvn et al.l 
( 1979allbl ) investigated steady flows across a range of temperature gradients and 



found that ionization equilibrium in the transition region is an acceptable assump- 
tion for iron at flow speeds no greater than 20 km/s, but that carbon and oxygen 
ion populations can be driv en aw ay from equilibrium a t flow speeds of only 1 km/s. 
iRavmond fc Dupred (|l978l) and IDup ree et all (|l979j) carried out a similar study 



related to steady flows in the transition region and also found sign ificant depar- 
tures of the ion populations from equilibrium. iBorini &: Noel ( 1982f ) investigated 



the ionization state in coronal loops in the temperature range 0.2 MK to 2 MK 
and showed that considerable deviations from equilibrium ionization could arise 
in average to low intensity loops characterised by high-speed fiows. They reported 
a pronounced effect for cooler loops, which despite exhibiting lower speed flows 
we re found to have s teeper temperature and density gradients than hot loops. 
iNoci et al] ( 1989f ) calculate the number density of carbon ions for a selection 



of coronal loops models in the case of steady-state, sub-sonic fiows (siphon fiows) 
and found departures from equilibrium of the ionization state for fiows of only a 
few km /s at the loop apex and for a factor of 10 slower at the base of the transition 



Spadaro et al.l ( 1990ah calculated the spectral line profiles of carbon ions 



region. 

formed in the transition region that are commonly used in spectroscopi c diagnostic 
studies. They used the number densities of the carbon ions found by iNoci et aD 
and found predominantly blue-shifted emission lines, which could not be 
reconciled with observations that show both up- and down-fiows in the transition 
region. The absence of re d-shifted emission was attributed to the assumption of 
spatially uniform heating. ISpadaro et al.l ( 1990bj) focused on the corona and cal- 



culated the emissivities of carbon and oxygen, both in and out of equilibrium, 
and found substantial differences between the resulting radiative loss functions. 
In the case of up-fiows (down-fiows) the radiative losses were generally enhanced 
(suppressed). One may understand this by considering an ion of relatively low 
charge state transported into a region of temperature significantly higher than the 
formation temperature of the ion in equilibrium; the ion will tend to emit more 
strongly since a greater proportion of the electrons in the bulk of the distribution 
wil l have sufficient energy to ex cite its emission lines. 

ISpad aro fc Ventural (|l994allbh studied the effect of non-equilibrium ion popu- 



lations on the line intensities of O VI and H I ions that originate in solar wind 
source regions. They calculated the intensity and line profiles for equilibrium and 
non-equilibrium ionization balance based on a steady fiow model, finding signif- 
icant deviations from equilibrium beyond 3 — 4 solar radii for O VI and beyond 
5 solar radii for the Lyman a emission from H I. These results are significant for 
estimates of the solar wind speed that rely on the Doppler-dimming technique, 
which estimates the speed from variations in the line intensities compared with 
the ir expected values in the absence of a steady outfiow. 

I Spadaro et al.l ( 19941) investigated the signatures that may be observable when 



non-equilibrium ion populations are present and should be considered when car- 
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rying out spectroscopic diagnostics using line ratios. Since non-equilibrium ion 
populations are displaced from their temperatures of peak abundance in equilib- 
rium, the temperature-dependent Boltzmann factors that appear in the expression 
for the excitation rate coefficients for the spectral lines are changed considerably, 
which results in changes to the energy level populations, the line intensities and, 
consequently, the values of the line ratios. The values of the temperature sensitive 
line intensiti es arising from th e non-equilib rium C IV and O IV- V I populations 
calculated bv lNoci et al.l ( 1989() (carbon) and lSpadaro et aT I (|l990bl) (oxy gen) were 
compared with the same line intensities computed in equilibrium. In the presence 
of a non-equilibrium ion population the line intensities were found to be reduced 
for both up-flows and down-flows across the transition region temperature gradi- 
ent. The C IV population was found to be the most sensitive to non-equilibrium 
ionization, with decreases in the line ratio by an ord er of magnitude in th e case of 
down- flows. In response to discrepancies identified bv lKeenan et al.l ( 19921) between 
C IV line intensities observed duri ng highly dynamic ev ents and theoretical predic- 
tions of the same line intensities. ISpadaro et afl (1995) used a siphon-flow model 
and non-equilibrium ion populations to recalculate the predicted line intensities. 
However, they found only a marginal improvements in the agreement between the 
observed and predicted intensities when non-equilibrium ionization was accounted 
for, and conlcuded that the observed int e nsitie s could not be reconciled with a 
sub-sonic, siphon-flow model. lEsser et~al ] (|199 8') examined the effect on the ion- 
ization state of the solar wind when the acceleration process occurs at much lower 
heights in the solar atmosphere than previously considered, based on flow speeds 
estimated from chromospheric, transition region and coronal emission lines. These 
observations yielded flow speeds for O VI ions that are a factor of 3 — 4 greater 
than indicated by earlier work, which imply the ion populations may depart from 
equilibrium as they are transported at speed across the steep temperature gradi- 
ents found in the lower atmosphere. In this case, the use of charge state ra tios to 
estimate equilibrium temperatures is unlikely to be valid. lEsser et al" ] ("1998) found 
outflow models with speeds in the regi on of 130 — 230 km/s t o predict charge state 
ratios consistent with those observed. lEdgar fc EsseH (|2000|) considered the effect 
of non-equilibrium ionization on the ratio of Ne VI to Mg VI lines in the solar 
transition region, which is used as a diagnostic of the first ionizaton potential 
(FIP) effect. In the presence of a strong heating or cooling effect the populations 
of ions of low FIP are enhanced relative to those of higher FIP. They calculated 
the non-equilibrium populations of these ions for simple fiows across the transi- 
tion region and showed that their spectral line ratios depend on non-equilibrium 
effects, as well as on the temperature and density. 

These investigations into the consequences of non-equilibrium ionization as- 
sumed steady-state conditions where only fiows may affect the ionization state in 
the presence of a steep temperature gradient. In general, this is due to the assump- 
tion of some form of constant heating that maintains the plasma in a steady-state 
condition, but this need not be so. There exist mechanisms by which energy can be 
impulsively released into the plasma on short timescales (e.g. a coUisionless shock 
or magnetic reconnection, (REFERENCES TO ARTICLES IN THIS VOL- 
UME TO BE ADDED)) leading to temperature changes on timescales that are 
short compared with the ionization time. Local temperature enhancements can 
give rise to localised pressure gradients which may in turn drive fiows. Conse- 
quently, a detailed understanding of the consequences of non-equilibrium ioniza- 
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tion requires a treatment of both local, temporal changes in the plasma properties 
and the fast transport of ions by flows. A shock is perhaps the simplest case, since 
it drives the plasma from one nearly steady state to a nother, and if th e shock is 
stong the ionization state can be far from equilibrium. iMa et alj ( 201l|) used the 
compression, density and heating determined from optical and radio observations 
of a CME-driven shock to compute the time-dependent ionization in the post- 
shock flow, and they found that it matched the observed rise times of emission in 
the AIA band s. 

iHansteenl (|l993l) presented a numerical model that solved the time-dependent 
ion population equations in tandem with the hydrodynamic equations, taking ac- 
count of departures from ionization equilibrium on the radiative losses for ions 
formed below 0.3 MK. The model was used to study the dynamic response of a 
coronal loop to energy released impulsively near the apex. It was found that the 
line shifts predicted for C IV, O IV and O VI by the model wer e consistent with th e 
persistent red-shifts observed in transition region lines (e.g. iBrekke et a"l]|l997l) . 
The amplitude of the predicted line shift was shown to depend on the ionization 
timescale of the emitting ion. It was also found that the radiative losses could 
change by a factor of 2 due t o the influence of flows and waves on the ion pop- 
ulation. I^riac^er^al] ( 1999a|) noted the presence of blue-shifts at temperatures 
characteristic of the transition region in the quiet Sun and in active regions and 
iTeriaca et al.' ('1999b') suggested that impulsive heating localised in the transition 
region at the temperature of peak O VI abundance in equiibrium (0.3 MK) might 
account for the presence of red-shifts and blue-shifts. The heating (whether located 
at the loop apex or in the transition region) generates compression waves and by 
including the partial reflection of the downward propagating wave from the chro- 
mosphere, and allowing for non-equilibrium ionization, reasonable agreement was 
found between the observed Doppler-shifts and those predicted by the numerical 
model (red-shifts i n the cooler C IV lines and blue-shifts in the warmer O VI lines) . 
iDovle et al.l (|2002l) found that the higher the temperature at which a heating event 
occurs then the greater the delay in the response from the mid-transition region 
lines in terms of changes in the Doppler-shift. 

Bradshaw & Mason (2003) studied the response of the plasma and the ioniza- 
tion state to a small-scale, impulsive energy release at the apex of a coronal loop, 
characteristic of nanoflare heating and solved the ion population equations for the 
15 most abundant elements of the solar atmosphere (including C, O, Ne, Mg, Si 
and Fe). The ionization state was used to calculate the radiative loss function in 
the energy equation, thereby coupling the energy balance with the ionization state. 
They concluded that broad / narrow-band imaging instruments can miss small- 
scale heating events entirely due to the weak sensitivity of the non-equiibrium 
emissivity to the changing temperature compared to the emissivity for equilib- 
rium ionization, which fell by a factor of up to 5. The non-equilibrium emission 
remained relatively steady throughout the heating event, despite a factor 2 change 
in the temperature on a timescale of 30 s. In order to diagnose non-equilibrium 
ionization they proposed searching for signatures in line ratios of ion pairs that 
are populous in the temperature range of interest but have different characteristic 
lifetimes (e.g. C IV and O VI in the transition region, or different ions of Fe at coro- 
nal temperature). Brad shaw et al 

] (11004) investigated non-e quilibrium ionization 
i n a s mall compact flare, using the same numerical model as lBradshaw &: MasonI 
(|2003h . and localised the energy release in the corona to drive the flare evolution 
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Fig. 6 Ion population fractions along a heated coronal loop calculated in equilibrium and 
non-equilibrium. The highly charged states associated with the hottest plasmas in equilibrium 
are never reached in the non-equilibrium calculation before the onset of cooling. Image credit: 
iReale fc Orlandol l l200a) . 



by thermal conduction. During the impulsive phase they found the emissivities of 
He I, He H and C IV in the transition region to be strongly enhanced above their 
expected equilibrium values, which was then followed by a significant reduction 
leading to an increase in the amount of chromospheric plasma ablated into the 
corona (less energy radiated in the transition region leaves more energy available 
to drive ablation). During the initial energy release the charge state of the coro- 
nal ions was seen to evolve substantially out-of-equilibrium with the increasing 
temperature and line ratio measurements would yield plasma temperatures that 
are much greater than the formation temperature of the emitting ion. During the 
gradual phase the emissivity at transition region temperatures was suppressed re- 
leative to equilibrium with reduced downflow velocities, since the enthalpy flux did 
not have to work as hard to power the transition region, and commensurately in- 
creased radiative cooling time-scales. The flare emission as it would be detected by 
TRACE in its 171 A and 195 A wavelength bands was computed and it was found 
that the filter ratio technique can give reasonably good estimates of the plasma 
temperature in quiescence. However, when the populations of Fe VIII, Fe IX, Fe X 
and Fe XII exhibited non-equilibrium effects the temperatures derived from filter 
ratio measurements were unreliable. 

iBradshaw fc Cargilll (|2006|) and iReale fc Orlandol (|2008l) considered strong or 
'explosive' heating, on short timescales to high temperatures, in an initially rar- 
efied coronal loop atmosphere, to determine the consequences for the evolution of 
the ionization state. They found extremely strong departures from ionization equi- 
librium and concluded that for sufficiently short heating events the charge states 
characteristic of the highest temperatures reached (10 — 30 MK) could never be 
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created before the onset of fast cooling by thermal conduction and coronal filling 
by chromospheric ablation (Figure[6|). In consequence, the emission measure peaks 
at temperatures significantly lower than the peak temperature of the plasma and 
forward modeling emission in the wavelength range of Hinode-EIS showed that no 
'hot' (e.g. > 10 MK) component of the plasma would be detected. Heating models 
that assume ionization equilibrium predict such a hot component, but no obser- 
vational evidence has yet been found and so non-equilibrium ionization presents 
one possibility to reconcile observations with current theory. In the future, obser- 
vations of the solar X-ray continuum could be used to confirm, or otherwise, the 
presence of a hot component to the emission. The X-ray continuum is due mostly 
to bremsstrahlung from H and He, which is not sensitive to non-equilibrium ion- 
ization effects. There is a significant contribution from radiative recombination 
in over-ionized plasmas and radia tive recombination continua (RRC) are seen 
in a few X-ray binaries (Cyg X -3: IPaerels et al .1120001) and supernova remnants 
(IC 443: lYamaguchi et al. I l2009l) . In under-ionized plasmas, such as may be cre- 
ated in the case of rapid heating in the solar corona, the RRC is weak compared to 
the bremsstrahlung. In any case the RRC scales as eyipihv /ksT) and so provided 
the edges of the most abundant elem ents are avoided, the continuu m shape can 
be used to diagnose the temperature. iBradshaw h KlimchukI (|201l[) conducted a 
more extensive survey of the parameter space of energy release magnitudes and 
timescales, and carried out forward modeling to predict the emission that would 
be detected in the passbands of the recently launched SDO-AIA, in order to per- 
form a more detailed evaluation of the potential for non-equilibrium ionization 
to explain the high temperature part of the emission measure. The study led to 
several conclusions: (1) Deviations from equilibrium were found to be greatest for 
short-duration nanofiares at low initial coronal densities. (2) Hot emission lines 
were the most affected and could be suppressed to the point of invisibility. (3) 
For many of the heating scenarios considered the emission detected in several of 
the SDO-AIA channels (131, 193, and 211 A) was predicted to be dominated by 
warm, overdense, cooling plasma. (4) It was found to be difficult to avoid creating 
coronal loops that emit strongly at 1.5 MK and in the range 2-6 MK, which are the 
most commonly observed kind, for a broad range of nanoflare scenarios; the mere 
abundance of such loops does not help to constrain the heating parameter space. 
(5) TheFeXV (284.16 A) emission predicted by most of the models was about 10 
times brighte r than the predicted Ca XVI I (192.82 A) emission, consistent with 
observations. iBradshaw fc KlimchukI ( 201ll) concluded that small-scale, impulsive 
heating that induces non-equilibrium ionization leads to predictions for observable 
quantities that are entirely consistent with what is actually observed. 

On larger spatial scales [Rakowski et al.l ( 20071) examined the ionization state 
of several elements derived from in-situ observations of a halo coronal mass ejec- 
tion (CME) made by the Advanced Composition Explorer (ACE). They assumed 
an evolution for the CME based on observations and models, and solved the ion 
population equations for the elements to be compared with the ACE data. They 
found that plasma in the core of the CME required further heating, possibly due 
to post-eruptive reconnection following the filament eruption, to reconcile the pre- 
dicted and observed populations. Plasma in the CME cavity, however, was found 
not to be further ionized following the eruption, because the low density in that 
region e ffectively freezes the io n po pulations in t he sta te they existed in close to 
the Sun. iMurphv et al.l ( 201l|) and iLandi et al.l ( 2010l) used the time-dependent 
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ionization in CME ejecta to constrain the temperature history of the expanding 
plasma and show that an amount of heat comparable to the kinetic en ergy must 
be inj ected to counteract the radiative and adiabatic expansion cooling. IKo" et al.l 
( 20101 ) examined the time-dependent ionization in post-eruption current sheets for 
a Petschek-type reconnection exhaust, and found that the observable line inten- 
sities depended strongly on the height of the reconnection X-line, as well as the 
density and magnetic field in the surrounding plasma. 



7.2 Non-Maxwellian Distributions 

The heat flux is sensitive to the underyling distribution of the particles that carry 
it and since properties of the distribution may also manifest in the emission spec- 
tra, then one may consider predicting observable signatures of the heat flux based 
on the results of numerical experimen ts and then searching for them in obser- 



vational d atasets. Karpen et al. (|l989|] used the results of their earlier flare cal- 



culations (jKarpen &: DeVordll987f) to investigate the effect of different heat flux 
formulations on the X-ray resonance lines of Ca XIX and Mg XI and included non- 
equilibrium ionization in their computation of the spectral lines. By comparing the 
results of spectroscopic diagnostics carried out with the predicted and observed 
emission lines, they found the properties of the flare plasma to be most consistent 
with the non-local formulation of the heat flux. Ijiang et al.l (^QQ§) found thermal 
conductivity supressed relative to the classical value in a loop-top source during the 
late decay phase of a flare, strongly indicative of the onset of flux limiting. The 
measured cooling timescale was longer than that predicted by classical thermal 
conduction, but shorter than for radiative cooling. However, they were unable to 
definitively determine whether plasma wave turbulence was providing additional 
he ating and / or surpres sing conduction by scattering electrons. 



lEsser fc Edeail (|200dl address the issue that electron temperatures observed at 
the solar wind acceleration site in the inner corona are too low to give rise to the ion 
populations observed in-situ in the solar wind, by considering non-Maxwellian elec- 
tron distributions. They show that reconciling the low electron temperatures and 
the relatively highly charged ions requires a number of conditions to be satisfied in 
the inner corona. (1) The electron distribution function must be near-Maxwellian 
at the coronal base. (2) A departure from Maxwellian must then occur rapidly as a 
function of height, reaching essentially interplanetary properties within a few solar 
radii. (3) Ions of different elements must have different speeds to separate their 
freezing-in distances enough that they encounter different distributions. They also 
show that the required distributions are very sensitive to the electron temperature, 
density, and ion fiow speed profiles in the coronal region where the ions form. 

A number of studies have adopted forms of non-Maxwellian electron distri- 
butions to ex plain differences bet ween the predicted and observed properties of 
emission lines. [Pufton et aD ( 19841) found that discrepancies between the observed 



and theoretically predicted ratios of lines from Si I II could plausibly be ex plained 
by non-Maxwellian electron distributions. However. [Anderson et all ( 1996j) showed 



that the k=2.5 distribution used in velocity filtration models of coronal heating 
over predicts the in t ensiti es of lines normally formed near lO'' K by a factor of 
100. iPinfield et all ( 19991) presented evidence that observations made by SoHO- 



SUMER showing enhancements in the predicted Si III 1313 A line intensity by a 
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factor of 5 in active regions, and by a factor o f 2 in the quiet Sun and coronal holes, 
could also be explained by non-Maxwellians. iRalchenko et aD ( 2007|) showed that 
the excess brightness of some hotter lines (low lying transitions in ions formed at 
temperatures greater than 2 MK) in the quiet corona may be accounted for by 
a two-component Maxwellian, where a high-energy component adding 5% elec- 
trons in the temperature range 300 — 400 eV is needed to account for the excess 
brightness. 



iDzifcakova et all (|2008cl) carried out thermal and non-thermal diagnostics of 
a solar flare observed with RESIK and RHESSI. They found that in compari- 
son with a synthetic isothermal or multithermal spectra, a non-thermal synthetic 
spectrum fitted the observed Si XII dielectronic satellite lines much more closely 
(with error less than 10%), and concluded that evidence for significant deviations 
of the free electron distribution from Maxwellian during the im pulsive phase of a 
solar fi are can be diagnosed using X-ray spectral observations. IDzifcakova et all 
( 2011a|) explained features of the RESIK X-ray fiare spectra using a Maxwellian 
or n-distribution for the bulk and a power-law tail, finding that the power-law 
tail has only a small effect on the satellite-to-allowed Si Xlld / Si XIII ratio, 
which is sensitive t o the shape of the bul k distribution and allows the parameter 
n to be diagnosed. iKulinova et all ( 201lh carried out diagnostics of non-thermal 
distributions in solar fiares observed with RESIK and RHESSI. They used two 
independent diagnostic methods, both indicating the fiare plasma affected by the 
electron beam can have a non-thermal component in the 25 keV range, which was 
found to be well-described by the n-distribution. Spectral line analysis also revealed 
that the n-distribution does not occupy th e same spatial locatio n as the thermal 
component detected by RHESSI at 10 keV. iKarlickv et al.l ( 20121) investigated the 
physical meaning of n-distributions in solar fiares. The electron component of the 
return current in a beam-plasma system was shown to have the form of a moving 
Maxwellian and this was found to be very similar to the high-energy part of an 
n-distribution. 



iDudfk et al.l (j2009l) calculated TRACE EUV filter responses to emission arising 
from non-Maxwellian distributions and showed that for K-distributions the result- 
ing responses to emission are more broadly dependent on temp erature, and their 

maxima are fiatter than for the Maxwellian electron distribution. IDzifcakova &; Kulinoval 
([2010) computed a set of synthetic spectra for various ^-distributions with varying 
electron densities and mean energies in the spectral range corresponding to the 
Hinode/EIS and Coronas-F/SPIRIT detectors. Strong EUV lines of Fe in various 
degrees of ionization were used to analyze the sensitivity of the line ratios to the 
shape of the distribution function, electron density, and temperature. It was found 
that EUV coronal Fe lines are generally not very suitable for diagnosing the non- 
thermal distributions due to their high sensitivity to electron density, but pairs of 
Fe XVII lines were reasonably good candidates for non-thermal diagnostics. Fi- 
nally, Dzifcakova & Kulinova (20 lib') was able to explain the observed intensity 
of the Si III spectrum in coronal holes, the quiet Sun and active region transition 
regions by adopting an underlying K-distribution for the electrons. 
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8 Summary and Future Directions 

We have reviewed a number of the microphysical processes occuring in optically- 
thin, astrophysical plasma environments, such as supernova remnants and the solar 
corona, that specifically influence their internal physics such as energy transport 
and atomic processes and, in consequence, their emission spectra. In particular, 
we have described the theory of spectral line formation in the coronal approxi- 
mation and how it is affected by the de-coupling of the ion population from the 
local temperature (non-equilibrium ionization), that arises when coUisional pro- 
cesses are unable to keep pace with heating or cooling, and by the formation of 
non-Maxwellian particle distribution functions, which are also related to the coUi- 
sionality of the plasma. Calculations of the charge state of the plasma both in and 
out of equilibrium, and the most computationally tractable formulations of the 
kinetic equations that are solved to find the underlying particle distribution func- 
tion have been presented. A selection of methods have been discussed by which 
the heat flux can be calculated in a hot plasma, when even near-thermal electrons 
have long mean-free-paths in relation to the characteristic spatial scales of the sys- 
tem, inducing strongly non-Maxwellian distributions, without recourse to solving 
a full kinetic equation. In addition, the ways in which non-Maxwellian distribu- 
tions alter the rates of excitation and ionization have been considered. The results 
and findings from practical applications of these calculations have been shown 
throughout and the physics of these processes has been connected with the total 
radiative energy loss from the system. Finally, evidence for observational signa- 
tures of non-equilibrium ionization and non-Maxwellian particle distributions has 
been presented in association with discussions of the diagnostics that have been 
used to reveal their infiuence. 

In the future, as astrophysical plasmas are probed with ever greater spatial, 
temporal and spectral resolution, we expect the microphysical processes that we 
have discussed here to become increasingly important to developing a full under- 
standing of the physics that drives and gover ns these systems. The forthcoming 
Interface Region Imaging Spectrograph (IRIS: lOe Pontieull2009[ ) will provide de- 
tailed spectroscopic observations of exactly the region of the Sun's atmosphere 
where non-equilibrium ionization and non-Maxwellian particle distributions may 
play a large role in forming the spectral emission lines and, consequently, in deter- 
mining what information about those regions can be extracted from the spectra by 
diagnostic studies. The upper-chromosphere and lower-transition region are highly 
dynamic environments where non-equilibrium ionization has been predicted to be 
a significant factor in emission from ions such as C IV and Si IV (Li-like and Na- 
like, respectively) that are undergoing heating ( Judge et al.ll201^ : their enhanced 



emission would lead to over-estimates of the density if equilibrium ionization were 
assumed when interpreting the observations. IRIS may be able to shed light on 
the strength of departures from equilibrium in this regime. 

Furthermore, streaming particles that enhance the tails of the particle distri- 
butions in the interface region may play a role in producing emission from ions of 
higher charge state than would be predicted from the local temperature alone. A 
larger tail population provides more electrons with sufficient energy to ionize the 
ambient plasma to a greater degree. The source of the streaming particles may be 
a hot (> 10 MK) component of the coronal emission due to in-situ heating in a 
high-altitude region where the energy per particle is large, leading to high tem- 
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peratures and mean-free-paths of lengths on the order of the spatial scales of the 
magnetic structures, even for near-thermal electrons. Observations of the corona at 
the highest spatial resolution so far achieved (75 km: ICirtain et al.|[2013[ ) indicate 
the presence of entwined bundles of magnetic flux that may be reminiscent of the 
long -theorized bra iding of magnetic field lines leading to reconnection and heating 
fe.g. |Parkeilll983h . The hot component may be extremely di fficult, if not impos- 
sible, to observe directly if the plasma is initially tenuous ( Bradshaw fc Cargilll 



120061: iReale fc OrlandolboOSi) . but if reconnection does lead to high temperatures 
and streaming particles, or direct particle acceleration, in the corona, then suit- 
ably sophisticated numerical models may be able to predict detectable signatures 
from lower altitude emission as indirect evidence that can be searched for in real 
observations by instruments such as IRIS. 

Finally, the initial state of what ultimately becomes coronal plasma after heat- 
ing occurs is another question that is worthy of attention. For example, does 
heating occur during active region emergence or following the draining of the ma- 
terial dredged up from below the surface as it rises? If the plasma carried to high 
altitudes cools below ^ 20,000 K then it becomes significantly partially ionized 
and, no longer supported by the Lorentz force due to the emerging field, the neu- 
tral atoms rain back onto the surface. In this scenario one might expect the active 
region plasma to be in an initially tenuous state and the energy per particle rel- 
atively high in the case of direct heating, giving rise to very high temperatures 
where the magnetic field strength and free energy are greatest (in the core of 
the active region). This may also be the case if the heating is intermittent and 
the corona is allowed to drain substantially between heating events. Evidence is 
beginning to accumulate to suggest that the frequency o f heating in active re- 
gions increases with its age ( Ugarte-Urra fc Warren 20121): voung active regions 



are h eated by low-frequency events (e.g. iMulu-Moore et al.ll201ll : lBradshaw et al 



2012[) and older active r egions are heated with greater frequency (e.g. lWarren et al 
201O; 'Ree p et al.ll20i;t ). The physics of heat fiux saturation and non-local ther- 
mal conduction must come into play when considering the energy transport and 
the overall energy balance of a hot but tenuous atmosphere and the treatments 
that extend the classical heat fiux, described in Section 15.21 must be revisited. 
In the case of steady heating, where the atmosphere is near hydrostatic, the fiux 
saturation regime is not reached but non-local thermal conduction may still be 
important in the high-temperature cores of active regions. 
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